function [ xdot ] = model(t, x)
%MODEL Summary of this function goes here
%   Detailed explanation goes here

global K productionRateFnc;

xdot = zeros(2, 1);
xdot(1) = productionRateFnc(t, x(1)) - x(2);
xdot(2) = -K * productionRateFnc(t, x(1));

end

